IFUP-TH 6/96 



A strong-coupling analysis of two-dimensional 0{N) a models 
with AT < 2 on square, triangular and honeycomb lattices 



Massimo Campostrini, Andrea Pelissetto, Paolo Rossi, and Ettore Vicari 
Dipartimento di Fisica deU'Universita and I.N.F.N., 1-56126 Pisa, Italy 



Abstract 

The critical behavior of two-dimensional O(A^) a models with < 2 on the 
square, triangular, and honeycomb lattices is investigated by an analysis of 
the strong-coupling expansion of the two-point fundamental Green's function 
G{x), calculated up to 21st order on the square lattice, 15th order on the 
triangular lattice, and 30th order on the honeycomb lattice. 

For N < 2 the critical behavior is of power-law type, and the exponents 
7 and v extracted from our strong-coupling analysis confirm exact results 
derived assuming universality with solvable solid-on-solid models. 

At = 2, i.e., for the 2-d XY model, the results from all lattices con- 
sidered are consistent with the Kosterlitz-Thouless exponential approach to 
criticality, characterized by an exponent o" = ^, and with universality. The 
value o" = I is confirmed within an uncertainty of few per cent. The prediction 
rj = J is also roughly verified. 

For various values of < 2, we determine some ratios of amplitudes 
concerning the two-point function G{x) in the critical limit of the symmetric 
phase. This analysis shows that the low-momentum behavior of G{x) in the 
critical region is essentially Gaussian at all values of < 2. New exact 
results for the long-distance behavior of G{x) when A^ = 1 (Ising model in 
the strong-coupling phase) confirm this statement. 

FAGS numbers: 75.10 Hk, 05.50.-hq, 11.10 Kk, 64.60 Fr. 
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I. INTRODUCTION 



The strong-coupling expansion is one of the most successful approaches to the study of 
critical phenomena. Many important results concerning physical models at criticality have 
been obtained by deducing the asymptotic critical behavior of physical quantities from their 
strong-coupling series. 

We have calculated the two-point Green's function 

G{x) = (4 ■ So), (1) 

of two-dimensional O(A^) a models on the square, triangular and honeycomb lattices, respec- 
tively up to 21st, 15th, and 30th order in the strong-coupling expansion. Such calculations 
were performed within the nearest-neighbor lattice formulation, described by the action 

5 = -iV/?^4,-4., (2) 

links 

where is a iV-component vector, the sum runs over all the links, and xi,Xr indicate 
the sites at the ends of each link. The comparison of results from strong-coupling series 
calculated on different lattices offers the possibility of important tests of universality, which, 
if positive, strongly confirm the reliability of the final results. 

A complete presentation of our strong-coupling computations for O(A^) a models in two 
and three dimensions will be presented in a forthcoming paper. A preliminary report on 
our calculations can be found in Ref. On the square lattice our strong-coupling series 
represent a considerable extension of the 14th order calculations of Ref. 0, performed by 
means of a linked cluster expansion, which have been rielaborated and analyzed in Ref. |Q . 
We also mention recent works where the linked cluster expansion technique has been further 
developed and calculations of series up to 18th order [Q and 19th order [^] for bulk quantities 
in (i = 2, 3, 4 have been announced. 

In this paper we focus on 2-d 0{N) a models with N < 2. The analysis of our strong- 
coupling series for models with > 3, i.e., those enjoying asymptotic freedom, is presented 
in Ref. ||. 

Two-dimensional O(A^) a models with N < 2 should present a standard power-law 
critical behavior, and should be described at criticality by conformal field theories with 
central charge c < 1. The most physically relevant models in this range of values of A^ are 
self-avoiding random walk models and Ising models, corresponding respectively to A^ = 
and A^ = 1. At A^ = —2 the critical theory has been proven to be Gaussian 0], i.e., 7 = 1, 
z/ = |, and 1] = 0. Assuming universality with solvable solid-on-solid models, exact formulas 



for the critical exponents in the range —2 < N < 2 have been proposed P|-pi][ , interpolating 
the critical behaviors at A^ = —2, 0, 1, 2. The critical exponents of the magnetic susceptibility 
7 and correlation length u are conjectured to be 

3 + a2 
^ ~ 4a(2 - a) ' 



2 



where the parameter a is determined by the equation 

'2n 



N 



-2 cos 



(4) 



with the constraint 1 < a < 2. The exponent rj can be obtained by the hyperscahng relation 
7 = (2 - r])iy. 

In the hmit N ^ 2 formulas (|^) yield 7 ^ oo and z/ — > 00, suggesting that at = 2 
the critical pattern should not follow a power-law behavior. The XY spin model in two 
dimensions, i.e., the N = 2 model, is conjectured to experience a Kosterlitz-Thouless phase 
transition [jll| of infinite order, characterized by a very weak singularity in the free energy 
and an exponential divergence of the correlation length at a finite f3. This model should 
describe the critical properties of a number of two-dimensional systems, such as thin films 
of superfluid helium. 

According to the Kosterlitz-Thouless (K-T) scenario, the correlation length is expected 
to behave like 



(5) 



for < r = 1 — P/ Pc ^ 1- The value of the exponent is cr = | and 6 is a non- universal 

00 of the 



positive constant. At the critical temperature, the asymptotic behavior for r 
two-point correlation function should be (cfr. e.g. Ref. fl^ ) 
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with = I and 9 = j^. Near criticality, i.e., 
susceptibility can be deduced from Eq. (10): 



for < r -C 1, the behavior of the magnetic 
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[l + 0(r'^lnr)]. (7) 



In addition, the 2-d XY model is characterized by a line of critical points, starting from 
(3 = (3c and extending to /5 = cxo, with r] going to zero as for (3 00. At criticality the 
2-d XY model should give rise to a conformal field theory with c = 1. 

Numerical studies based on Monte Carlo simulation techniques and high-temperature 
expansions seem to support the K-T behavior, but a direct accurate verification of all the 
K-T predictions is still missing. As pointed out in Ref. |T3|, for j3 < Pc the corrections to 
the asymptotic behavior (§) should become really negligible only at very large correlation 
lengths, out of the reach of standard Monte Carlo simulations on today's computers, which 
allow ^ ^ 100 (cfr. Ref. |T^, where simulations for correlation lengths up to ^ ~ 70 were 
performed on lattices up to 512^). Monte Carlo simulations supplemented with finite-size 
scaling techniques allow to obtain data for larger ^. Ref. |lH] shows data up to ^ ^ 850, 



which, although consistent with the K-T prediction, do not really exclude a standard power- 
law behavior.]^ 



^Actually the author of Ref. [15| claims to favor a conventional power behavior to explain some 
discrepancies in the determination of the critical exponent 77. 
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Finite-size scaling investigations at criticality are required to be very precise in order to 
pinpoint the logarithm in the two-point Green's function. On the other hand, if this logarith- 
mic correction is neglected, the precise check of the prediction t] = ^ at (3c may be quite hard. 
The relevance of such logarithmic corrections and some of the consequences of neglecting 
them have been examined in Ref. [0. Numerical studies by Monte Carlo renormalization- 
group and finite-size scaling techniques seem to favor a lower value of 1], which might 

be caused by the neglected logarithm. The most accurate verification of the K-T critical 
pattern has been shown in Ref. |jl3| by numerically matching the renormalization-group 
trajectory of the dual of the XY model with that of a the body-centered solid-on-solid 
model, which has been proven to exhibit a K-T transition. The advantage of this strategy 
is that such a matching occurs much earlier than the onset of the asymptotic regime, where 
numerical simulations can provide quite accurate results. 



The analysis of the strong-coupling series (cfr. e.g. Refs. |1T8|JI9|1 , where a few moments 



of the two-point Green's function were calculated on the square and triangular lattices 
respectively up to 20th and 14th order for the special value N = 2) substantially supports 
the K-T mechanism, but it does not provide precise estimates for the exponents a, rj, and 
6, probably for two reasons: (i) the asymptotic regime in the terms of the series may be set 
at very large orders; (ii) the logarithmic correction may cause systematic errors in most of 
the analysis employed. 

The computation of strong-coupling series on the honeycomb lattice and the extension of 
series on the square and triangular lattices motivate a new strong-coupling analysis of the 2-d 
XY model. We focus on the K-T mechanism, searching for evidences of this phenomenon. 

As already shown in Refs. P,pO|, the strong-coupling analysis provides quite accurate 
continuum-limit estimates when applied to dimensionless ratios of universal quantities, even 
in the case of asymptotically free models, i.e., when the critical point is = oo- We define 
some dimensionless ratios of scaling quantities (ratios of amplitudes) which characterize the 
low momentum behavior of the two-point function G{x), and estimate their values in the 
critical regime by directly analyzing their strong-coupling series. This will allow us to check 
how close the low momentum critical behavior of G{x) is to Gaussian behavior. 

The paper is organized as follows: 

In Sec. m we investigate the critical behavior of 2-d 0{N) a models with < 2 on 
the square, triangular and honeycomb lattices, extracting the relevant critical parameters 
by the analysis of the available strong-coupling series. For < 2 we compare the strong- 
coupling estimates of the critical exponents with Eqs. (^. For N = 2, i.e., the 2-d XY 
model, we verify the predictions of the K-T critical theory. In particular Sec. |II A| presents 
the general features of our strong-coupling analysis. Sees. |11 B| , p.1 L] , [11 D| and [11 F| contain 
detailed reports of the derivations of the various results; they are rather technical and can 
be skipped by readers not interested in the details of the analysis. In Sec. [II F[ the principal 
results are summarized and some conclusions are drawn. 

In Sec. m we evaluate, at criticality, the values of some amplitude ratios concerning the 
low-momentum behavior of G{x). We will present results for the most physically relevant 
models with N < 2, i.e., those with A^ = 0, 1,2. We also discuss their implications on the 
low-momentum behavior of G{x) in the critical region of the symmetric phase. 

In Sec. 1^ some new exact results concerning the asymptotic large-distance behavior of 
G{x) for the Ising models on the square, triangular and honeycomb lattices are presented. 
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In Apps. ^ |B], and we present, for = 0,1,2, the strong-coupling series of some 
relevant quantities used in this study, respectively for the square, triangular, and honeycomb 
lattice. 



II. STRONG-COUPLING ANALYSIS OF THE CRITICAL BEHAVIOR 

A. Analysis of the series 

From the Green's function G{x) one can derive many interesting quantities. Defining 
the moments of G{x) 

m,^=J2{xyG{x), (8) 

X 

we computed on each lattice the magnetic susceptibility and the second moment correla- 
tion length ^Q, 

^2 = ^ (9) 

Models with N < 2 should have a power-law critical behavior, which may be appropri- 
ately investigated by analyzing the strong-coupling series of x and in order to extract 
the critical exponents 7 and u. For = 2, in order to check the exponential approach 
to criticality predicted by the K-T mechanism and extract the relevant exponent a, as in 
Ref. |jl8[, we consider the strong-coupling series of the logarithm of x and ^g- More precisely, 
since x = 1 + 0{P) and = \c(3 + 0{P'^), where c is the coordination number of the lattice 
(c = 4,6,3 respectively for the square, triangular and honeycomb lattice), we consider the 
series 

00 

/;, = rMnx = Ec#, 

l^^p-Hn(^]=f:d.p\ (10) 



i=0 



According to Eqs. (|^) and (|^ and should behave as 

k-k- ^-"^ (11) 

and are therefore suitable for a standard analysis by Fade or integral approximants. A 
vanishing exponent a would indicate a standard power-law critical behavior. Conversely a 
stable non-zero value of a would exclude a power-law behavior. 

Estimates of the critical exponents can be obtained by employing the so-called critical 
point renormalization method (CFRM) The idea is that, when 

i 

B{x) = J2hx'--{xo-x)-^, (12) 
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we have 

h 



= ^ -x^ ~ (1 - (13) 

where now the position of the singularity is known. Therefore the analysis of the series C 
may provide an unbiased estimate for the difference between the critical exponents of the two 
functions A and B. In particular this idea can be applied to the case B = A^, allowing one 
to get a direct estimate of the critical exponent of A, provided a sufficiently large number of 
terms is known. The reliability of the determination of the critical exponent by this method 
may be checked by comparing the results for the critical point of an unbiased analysis with 
the exact result Xc = I- 

A general technique to extract physical information from a nth order strong-coupling 
series S{x) = J27=o ^i^^ is constructing approximants A{x) such that 

A{x) - S{x) = 0(x"+^), (14) 



and studying their singularities. For a review on the resummation techniques cfr. Ref. p2[ 



[l/m] Fade approximants (FA's) are ratios of two polynomials of degree / and m respectively, 
such that their Taylor expansion is equal to S{x) up to 0{x^^™). FA's are expected to 
converge well to meromorphic analytic functions. More flexibility is achieved by constructing 
FA's of the logarithmic derivative of S{x) (Dlog-FA analysis), and therefore enlarging the 
class of functions which can be reproduced to those having singularities of the form (z — zo)'^- 
[l/m] Dlog-FA's are obtained by integrating the [l/m] FA's of the logarithmic derivative of 
S{x). Then a [l/m] FA uses n = I + m terms of the series, while a [l/m] Dlog-FA requires 
n = I + m + 1 terms. 



Other kind of approximants can be constructed as solutions of differential equations 
We consider integral approximants (lA's) obtained from a first order linear differential equa- 
tion 

QMfix) + Pi{x)f{x) + R,{x) = O (x'^+'+^+^X ^ (^5) 



where Qm{x), Pi{x) and Rk{x) are polynomials of order m, /, and k respectively, and we fix 
<5m(0) = 1. These approximants are singular at the zeroes Xq of Qm{.x), and behave as 

A{x)\x-Xfi\-^ + B{x), (16) 

where A{x) and B{x) are regular in the neighborhood of xq, and 

= (17) 

When we analyze a nth order series, m, / and k must satisfy the condition k + l + m + 2 < n . 
If the position of the singularity xq is known, such an analysis can be easily modified forcing 
the approximant to have a singularity at Xq by substituting Qm{,x) — (1 — x/xo)Qm(a;), 
where Qm{x) is still a polynomial of order m with (5m(0) = 1. 

Unlike Dlog-FA's, lA's are suited to take into account subdominant terms in the vicinity 
of singularities, thus reducing possible systematic errors in the resummation of the series. 



6 



On the other hand, in order to get stable and therefore acceptable results, lA's require in 
general more terms in the series to be resummed than PA's or Dlog-PA's. 

As a final estimate from each analysis we took the average of the results from quasi- 
diagonal (non-defective) approximants (PA's or lA's) using all available terms of the series. 
The errors we display are just indicative, and should give an idea of the spread of the results 
coming from the various approximants which can be constructed from the series at hand. 
They are the square root of the variance around the estimate of the results coming also from 
quasi- diagonal approximants constructed from shorter series by one and two terms. In the 
following we will specify the approximants considered in each analysis. This procedure does 
not always provide a reliable estimate of the systematic error, which may be underestimated 
especially when the structure of the function cannot be well reproduced by the class of 
approximants used. A more reliable estimate of the true uncertainty should come from the 
comparison of results of different analysis of the same series, and from the analyses of series 
of different estimators of the same quantity, which in general are not expected to have the 
same analytic structure. 



B. Critical behavior of models with —2 < AT < 2 

In order to determine the critical exponents 7 and v of 2-d 0{N) a models with N < 2, 
we analyze the strong-coupling series of x and on square (21st order), triangular (15th 
order) and honeycomb (30th order) lattices. For such models, an analysis of the 14**^ order 
strong-coupling series on the square lattice, calculated in Ref. has been done in Ref. 0. 

For = (the self- avoiding walk), longer series are available To compare with 

the literature, observe that we have pT 



G{x) = J2P^ci{x), (18) 



I 



where q(x) is the number of self-avoiding walks of length / going from to x. Therefore 

X = J2i(3^ci, and 

xeG = iT.p'c,{Ri)r (19) 

I 

where q = J2xCi{x) is the total number of self-avoiding walk of length / starting from the 
origin, and 

Rl)=-T.^'c,{x) (20) 

is called the "mean end-to-end distance" in self-avoiding walks literature. 

Table | shows the results of a Dlog-PA analysis, reporting, for various values of and 
for each lattice, P^^^ and 7 as obtained from the strong-coupling series of x, and P^^^ and u 
from that of ^q. Differences between P^^^ and should give an idea of the real uncertainty 
on Pc- In the analysis of x we considered Dlog-PA's with I + m > 18 and m > / > 8 on 
the square lattice, / + m > 12 and m > / > 5 on the triangular lattice, / + m > 27 and 
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m > / > 12 on the honeycomb lattice. In the analysis of P'^Cg considered Dlog-PA's 
with I + m > 17 and m > / > 8 on the square lattice, I + m > 11 and m > / > 5 on the 
triangular lattice, Z + m > 26 and m > / > 12 on the honeycomb lattice. We tried also lA's, 
obtaining consistent results, which however only in few cases turned out to be more precise 
than those of the Dlog-PA's, so we do not report them. For sake of completeness and also 
to give an idea of the precision we can achieve with such an analysis, in Table | we report 
also results for = 0, 1 as obtained from our series, although exact results independent of 
the conjecture (0) exist for such values of A^. We warn that the errors displayed in Table | 
are related to the spread of the results from the Dlog-PA's considered, according to the 
procedure described in the Sec. p.1 A| , and therefore they are not always reliable estimates of 
the uncertainty. 

In the range — 1 ^ A^ ^ | the formulas (|^) for the exponents 7 and u are well reproduced, 
and universality is verified. Less precise determinations are obtained when approaching the 
endpoints A^ = ±2, presumably due to a rather slow convergence of the corresponding series 
to their asymptotic regime. 

We note that for models with A^ ^ 1 on the honeycomb lattice the physical critical 
point is not the singularity closest to the origin, but there is a pair of closer singularities 
on the imaginary axis. For instance, in the Ising model the physical singularity is placed 
at /3c = 1(2 + 73) = 0.658478... and there is a pair of singularities at /? = ±z7r/6 [p7[ 



Nevertheless Dlog-PA's of the magnetic susceptibility reproduce the physical singularity very 
precisely: the [15/15] Dlog-PA gives Pc = 0.658480 and 7 = 1.74993, to be compared with 
the exact result 7 = |- The unphysical singularities can be mapped away from the origin by 
performing the change of variable /3 — > z = tanh/3, where z is the character coefficient of the 
fundamental representation. With decreasing A^, /3c decreases while the above-mentioned 
imaginary singularities move away, so that at A^ ^ 0, the singularity closest to the origin is 
on the real axis, i.e., it is the physical critical point. 

For later comparison with the strong-coupling analysis of 2-d XY model, we have also 
analyzed the series of l^, defined in Eq. (0), for the Ising model on the square lattice. Since 
the critical behavior is of power-law type, should have a logarithmic singularity at /3c, and 



therefore an analysis like lA should give a ~ p3[. Indeed most of the lA's of the 20th order 
series of give \a\ ^ 0.02. We mention that a Dlog-PA analysis leads to misleading results 
in this case, since in order to reproduce the logarithmic behavior it gives rise to spurious 
singularities, which make the estimate of a at (3c unreliable. 



C. The 2-d XY model on the square lattice 

On the square lattice a strong-coupling analysis of the lowest moments of G{x) evaluated 
up to 20th order can be found in Ref. |1^. Having achieved further extension by one 



term of such series, we update here the situation on the square lattice. We note that the 
series obtained from our calculations (some of them are reported in App. ^ present small 
discrepancies with those reported in Ref. ITH]: they are in agreement up to 16th order. 
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but slight different at higher ordeiQ. We are confident that our series are exact, since they 
were generated for arbitrary and we have checked their N oo hmit against the exact 
solution, and compared them with the existing series for = 0, 1. 

We analyzed the 20th order series of = j3~^ Inx by both Dlog-PA's and lA's. We 
found (3c = 0.560(2) and a = 0.53(4) from Dlog-PA's (with / + m > 18 and m > / > 7). The 
integral approximant analysis, whose details are given in Table |I|, leads to jSc = 0.558(2) 
and a = 0.49(8) (considering lA's with m + I + k > 19 and m > l,k > 5). 

From the strong-coupling series of and we have constructed a new series ac- 
cording to CPRM, and analyzed it by standard methods: Dlog-PA's and lA's. We obtained 
a = 0.51(4) by Dlog-PA's (with I + m > 18 and m > I > 8) biased by imposing the pres- 
ence of a singularity at Xc = 1, and a = 0.50(2) by biased lA's (with m + I + k > 16 and 
m > l,k > 5). Table ^ shows some details of the lA analysis of the series of A^^-. From 
unbiased Dlog-PA's and lA's Xc is found to be equal to one within a few per mille, assuring 
us on the reliability of the estimates of the exponent a by this method. 

The above unbiased analyses strongly support the K-T prediction (|^). Although the 
estimate of a does not yet reach the high level of precision which is usually found in the 
analysis of strong-coupling series of considerable length, we can safely conclude that the 
value 0" = I is well verified with an uncertainty of less than 10% on the square lattice. 

Unbiased approximants of give less stable but definitely consistent results: we found 
a = 0.59(6) from CPRM, i.e., from an lA analysis biased at Xc = 1 (with 171 + 1 + k > 15 and 
m,l,k > 5) of the series Ag constructed from and /| according to CPRM. On the other 
hand by applying CPRM to the two series and one may verify that = a^. From an 
lA analysis biasing Xc = 1 (with m + l + k > 15 and m,l,k > 5) we found a^ — a^ = 0.010(6), 



which represents a good check of Eq. (11 



Once (y = \ was reasonably verified, we performed a set of biased analysis fixing a = | in 
order to determine the critical point. A way to bias the value of the exponent at cr = | is to 
analyze the series of the square of and by PA's. By doing so we obtained /5c = 0.5579(3) 
from and (5^ = 0.558(1) from Z| (we used PA's with / + m > 17 and m > I > 8). Still 
biasing a = |, lA's yield Pc = 0.5583(2) form l^, and Pc = 0.559(1) from (here we 
determine, in a lA analysis biasing the position of the singularity, the value of /3c which 
produces the exponent cr = |.) 

No complex singularities closer to the origin than /3c are detected in the various strong- 
coupling analysis, thus indicating that /3c is also the convergence radius of the strong-coupling 
expansion. 

For /? < /?c we have compared x and as obtained from our strong-coupling analysis 
with numerical data, available in the literature up to /5 = 1/1.96 = 0.5102... (corresponding 
to a correlation length ^ ~ 70) by using standard Monte Carlo simulations JI^, and up to 
P = 1/1.87 = 0.5347... by employing also finite-size scaling techniques |jl5|]. Actually most 
of the Monte Carlo data of ^ reported in Ref. [0] concern ^exp, i-e., the correlation length 
extracted from the long distance exponential behavior of G{x), but as we shall see in the 
next section ^c/^cxp — 0.999 at criticality. In order to get strong-coupling curves of x and 



^The difference is however small, at most 10 ^, and it does not change the conclusion of Ref. |18|. 
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as functions of (3, we used approximants constructed biasing cr = |, i.e., PA's of and /|, 
and lA's of and biasing cr = i. Fig. |l] compares strong-coupling curves of Inx with the 
available Monte Carlo data. Table [V| reports estimates of by strong-coupling expansion 
and Monte Carlo simulations for various values of /3. The agreement among the different 
calculations is satisfactory. 

In order to determine the exponent rj without biasing a, we considered various quantities 
which can be constructed using the lowest moments of the two-point Green's function: 

A = ^ = ^(l-^) + 0(r-lnr), (21) 



(22) 



C, = (/5/,)-Mn(^^j =| + 0(r-lnr). (23) 

An estimate of t] may then be obtained by resumming the corresponding strong-coupling 
series by PA's and Dlog-PA's and evaluating them at jSc- Since all of the above quantities 
are equally good estimators of the exponent t], differences in the results of their analysis 
should give an idea of the systematic error in the procedure. 

PA's and Dlog-PA's of (with Z -|- m > 16 and m > I > 7) lead to quite stable results: 
T] = 0.228(2) (where the error displayed, beside the spread of different approximants, takes 
into account the uncertainty on /3c). Similarly we found r] = 0.270(5) and r] = 0.226(5) 
respectively from PA's and Dlog-PA's of By (with I + m > 17 and m > I > 8) and of 
f3~^Crj (with Z -|- m > 15 and m > I > 7). The differences in such determinations indicate 
a systematic error of about 10%, and within 10% all results are consistent with the K-T 
prediction rj = ^. Furthermore, when analyzing the energy-series of (by performing the 
change of variable f3 ^ E and evaluating the corresponding approximants at Ec — 0.722, 



i.e., the energy value at Pc — 0.559 [Q), we obtained again a rather stable result but 
1] = 0.207(5), confirming the presence of a systematic error of about 10%. 

A source of systematic error in this analysis is the 0(r°" In r) correction expected in 
Eqs. ([2T|), (|2^) and (^31) which cannot be reproduced by PA's and Dlog-PA's. In particular 



Eq. (^Tj) implies a behavior 

DlogA^ ~ t"-^ (24) 

in the proximity of to jSc- In the Dlog-PA's the above singularities should be mimicked by a 
pole shifted at a /3 larger than /3c. Indeed in the analysis of the series of Ajj we have found 
a singularity typically at /5 ^ 1.1 1.2/5c. This fact will eventually affect the determination 
of Ar^ close to by a systematic error. However, since the singularity is integrable, the 
error must be finite, and the analysis shows that such errors are actually reasonably small. 
The behavior (pT]), neglecting the logarithm, could be reproduced by lA's , but we did not 
obtain sufficiently stable and therefore acceptable results by them. As a further check we 
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also performed a biased Dlog-PA analysis of the series of the quantity — |(1 — r^) ~ r*^ 
(neglecting logarithmic corrections). By fixing t] = ^ and (3c = 0.559 we found an exponent 
a ~ 0.6, which is satisfactorily close to the expected value cr = ^. 

The exponent 9 defined in Eqs. (|) and (|^ may be extracted by the analysis of the series 
of the ratios x/^g~^ and ^^4/^^"^, indeed 

^-^-r-^'-'il + Oir'^lnr)]. (25) 

Fixing = |, we performed biased analyses of the 20th order strong coupling series of the 
above ratios imposing /3c = 0.559. Dlog-PA's (with I + m > 17 and m,l > 8) provide 
the following estimates for 6 (obtained taking cr = |): 6 = —0.042(5) from x/^G~^y ciiid 
6 = —0.05(2) from 77^4/^^"'', where errors take into account, beside the spread of the Dlog- 
PA results, also the uncertainty on jSc- These numbers, although they confirm the fact that 
1^1 is small, are rather different from the K-T prediction 6 = j^. As already mentioned 
above, a source of systematic error for a Dlog-PA analysis is the correction 0(r^/^ In r) to 



the leading r~ behavior in formula (^Sf) . Since 6 is very small, this could cause a relevant 
departure from its true value. The analysis by lA's biased at Pc — 0.559 of the strong 
coupling series of the ratio x/^cf^ does not provide sufficiently stable results for 6. 



D. The 2-d XY model on the triangular lattice 

On the triangular lattice strong-coupling series of some lowest moments of G{x) have been 
calculated in Ref. up to 14th order. We calculated G{x) up to 15th order, thus extending 



by one order earlier calculations. We must again mention the existence of discrepancies 
between our calculations (cfr. App. P) and those of Ref. [19| in the 14th order terms (again 
of the order of 10~^). 

We performed on the triangular lattice the kind of analysis presented in the previous 
Subsection for the square lattice. We analyzed the 14th order series of = (3~^ Inx by both 
Dlog-PA's and lA's. We found (3^ = 0.3413(3), a = 0.52(2) from Dlog-PA's (with / + m > 11 
and m > / > 5), and (3^ = 0.33986(4), a = 0.473(3) from lA's with m + I + k > 12 
and m > l,k > 3 (the apparent stability of these lA determinations should not be taken 
seriously, we remind again that what it is really important for estimating the uncertainty is 
the comparison of results from different analysis). 

The CPRM applied to the strong-coupling series of and gives a = 0.53(1) by biased 
Dlog-PA's having a singularity fixed at Xc = 1 (with I + m > 12 and m > / > 5) and 
cr = 0.50(2) by biased lA's (with m + I + k > 11 and m > l,k > 3). Table |^ shows some 
details of the lA analysis. From unbiased Dlog-PA's and lA's, Xc is found to be equal to one 
within a few per mille, assuring us on the reliability of the estimates of the exponent a by 
this method. 

The critical point renormalization method gives good results also when applied to the 
series of leading to a = 0.52(4). By applying CPRM to the two series and one 
finds — (T^ = 0.001(2) from an lA analysis biased at Xc = 1 (with m + I + k > 12 and 
m > l,k > 3). 
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We can conclude that the K-T prediction o" = | is strongly supported by the above 
results. 

We also performed a set of biased analysis fixing a = | in order to determine the critical 
point. By analyzing and /| by PA's, we obtained (3c = 0.3400(2) and [3c = 0.3393(1) 
respectively (from PA's with I + m > 12 and m > / > 5). lA's of biased at a = | yield 
13c = 0.3410(5), while when applied to such analysis does not lead to relevant results, 
since it gives rise to spurious singularities in the real axis. A o" = | biased estimate of the 
critical point is then Pc = 0.340(1). 

We obtained estimates of t] by resumming the series of and i?^ (cfr. Eqs. ( [2T| ) and 
(0)) by PA's and Dlog-PA's and evaluating them at (3c ^ 0.340. PA's and Dlog-PA's of 
(with I + m > 10 and m > / > 5) and of 5^ (with / + m > 11 and m > / > 5) lead again 
to quite stable but shghtly discrepant results: respectively r] = 0.221(2) and r] = 0.270(4). 
The causes of possible systematic errors in the determination of r] are the same as for the 
square lattice, and we refer to Subs. |11 C| for a discussion. 

The analysis of the 14th order strong coupling series of x/^g"'' ~ ''"^^ biased by (3c — 0.340 
(using Dlog-PA's with / + m > 11 and Z,m > 5) yields the estimate 6 = —0.045(3), which 
is consistent with the square lattice result, but not with the K-T prediction. lA's does not 
provide sufficiently stable results also in this case. 

E. The 2-d XY model on the honeycomb lattice 

On the honeycomb lattice we calculated series longer than on the square lattice, up to 
30th order. Here the possibility of reaching larger orders is related to the smaller coordination 
number. However longer series do not necessarily mean that more precise results can be 
obtained from their analysis. This possibility is related to the approach to the asymptotic 
regime of the series, which is expected to be set later on lattice with smaller coordination 
number. Actually, as we shall see, the 30th order series on the honeycomb lattice provide 
results consistent with the K-T theory and universality, but less precise that those obtained 
from the series on the square and triangular lattices. 

Unbiased analyses of the series for lead to: (3c = 0.884(1) and cr = 0.55(1) from 
Dlog-PA's (with / -f m > 27 and m > / > 12), and (3c = 0.877(6), a = 0.4(2) from lA's 
(with m + I + k > 27 and m > l,k > 8). The stability of Dlog-PA's is suspect in this 
case, indeed we found that, just by adding to the series a simple constant of the order of 
unity, the change in the estimate of a turns out to be much larger than the error evaluated 
from the spread of the estimates of different approximants of the same series. Unlike PA's 
and lA's, the critical parameters provided by Dlog-PA's do not remain strictly invariant by 
adding a simple constant; this is only an asymptotic property of the Dlog-PA analysis. An 
analysis based on CPRM method fails to give stable results in this cases, probably because 
the available series are not sufficiently long to have their asymptotic regime set. 

More stable results are obtained when biasing the exponent at a = |. A check of 
consistency would be requiring that the critical points as extracted from and are the 
same. We obtained (3c = 0.879(1) from PA's (with / + m > 26 and m > / > 13) of and 
(3c = 0.878(2) from PA's of Z|, which is satisfactory. 

It is worth noticing that, unlike what happens on the square and triangular lattices. 
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on the honeycomb lattice the real singularity corresponding to the critical behavior of the 
theory is not the singularity closest to the origin in the complex /?-plane. A pair of imaginary 
singularities at /5 ~ ±i0.482 is detected in the analysis of the strong-coupling series of x- 
Taking as an estimate of the critical point (3c — 0.880, we evaluated rj from the series of 
and -B^ defined in Eqs. ( pT]) and (p2D. Again the value obtained from is about 10% 
lower than \'- Tj = 0.231(3) (from PA's and Dlog-PA's with / + m > 26 and m>l> 11), and 
that from is about 10% higher: rj = 0.28(1) (from PA's and Dlog-PA's with I + m > 27 
and m > I > 12). The behavior of the estimates from and i?^ observed in the various 
lattice seems to indicate that the source of systematic error is in a sense universal, i.e., it 
essentially depends on the quantity considered and approximately independent of the lattice. 

We again estimated the exponent 6 from the 29th order strong coupling series of the 
x/^g"''- Dlog-PA's (with l + m > 26 and /,m > 13) biased at /3c ~ 0.880 give 9 = -0.042(6), 
which is consistent with the estimates from the other lattices. 



F. Conclusions 

We have studied the critical properties of 2-d 0{N) a models with iV < 2 on the square, 
triangular and honeycomb lattices, by analyzing the strong-coupling expansion of the lowest 
moments of the two-point fundamental Green's function. 

The analysis of the strong-coupling series of x ^-iid ^ on the square, triangular and 
honeycomb lattices has substantially confirmed that models with N <2 present a standard 
power-law critical behavior with critical exponents given by Eqs. (^). We obtained rather 
precise determinations of the critical exponents in the region — 1 ^ ^ | (cfr. Table |), 
where formulas (||) are verified within 1%. The strong-coupling analysis becomes less precise 
approaching the endpoints N = ±2, presumably due to a rather slow convergence of the 
corresponding series to their asymptotic regime. Universality among models on the square, 
triangular and honeycomb lattices has been verified. 

The determinations of /5c and a for the 2-d XY model on the three different lattices are 
summarized in Table 0. These results are consistent with the K-T exponential approach to 
criticality and with universality. The best estimates of a come from the analysis of the series 
of the magnetic susceptibility, leading to a confirm of the value c" = | within an uncertainty 
of few per cent. The analysis of the series of the correlation length yields consistent 
results. The critical point renormalization method provides the most precise unbiased 
estimates of cr on the square and triangular lattices. These results rule out the possibility 
of a standard power-law critical behavior. 

On the square lattice most estimates of Pc yielded by our analyses he in the range 
0.558 < /5c < 0.560, although the lowest value Pc = 0.558 seems to be favored. This value is 



consistent with the results of an exponential fit to data of ^ up to ^ ~ 850 [|15l, which yielded 
Pc = 0.5593(13) and cr = 0.46(3), and with a biased exponential fit fixing cr = | to data up to 
^ ~ 70 produced by a standard Monte Carlo simulations [Q, which gave Pc = 0.559(3). But 
it is slightly smaller than the quite precise Monte Carlo renormalization group determination 
of Ref. [T^: Pc = 0.55985(25). The comparison of (suitable) resummations of the strong- 
coupling series of x and with Monte Carlo data (available up to ,^ ^ 850) turns out to be 
quite satisfactory (cfr. Table 0), giving further support to our conclusions. 
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The prediction i] = j is also substantially verified. By using different estimators of 77 we 
control the systematic error of our analysis, which turns out to be about 10%, and within 
about 10% our estimates of t] are always consistent with the value V = Z' 

Substantial discrepancies from the Kosterlitz-Thouless theory are found in the estimates 
of 6 we obtained on all lattices considered and using also different estimators. Our strong 
coupling analysis based on Dlog-PA's leads, similarly to the K-T prediction, to a small 
absolute value of 6, but it would favor the value 6 ~ —0.04, against the K-T value 6 = 
j^. Our strong coupling estimate seems to pass the universality check by changing lattice 
and estimator. On the other hand, we suggest some caution in considering our strong 
coupling estimate of 6. Given the smallness of its value, we cannot exclude that the observed 
discrepancy is due to systematic errors caused by the fact that Dlog-PA's cannot reproduce 
the correction 0(r^/^ In r) to the leading behavior in formula (|25|). Since this correction 
is expected to be present in all quantities we considered to estimate 6 (even when defined on 
different lattices), if its coefficients in the various cases are quantitatively similar the error 
might be about the same, and explain the apparent universality of our results. Moreover, 
the more general analysis based on lA's does not provide sufficiently stable results when 
applied to estimate 6, likely because the available series are not sufficiently long to this 
purpose. We mention that in Ref. an analysis based on Monte Carlo simulations led to 
the estimate 6 ~ 0.02, which is not consistent with both the K-T prediction and our strong 
coupling estimate. 



III. LOW-MOMENTUM BEHAVIOR OF G{x) IN THE CRITICAL REGION 

In this section we study the low-momentum behavior of the two-point fundamental 
Green's function in the critical limit of the symmetric phase. To this purpose, we consider 
the dimensionless renormalization-group invariant function 

Igl. 

In the critical region of the symmetric phase L{p, (3) is a function of the ratio y = p^/Mq 
only, where Mq = I/'Cg and C,g = '"^2/4x is the second moment correlation length, already 
introduced in the previous section. L{y) can be expanded in powers of y around y = 0: 

L{y) = l+y + l{y) 

00 

l{y) = Y^c.y\ (27) 

i=2 

l{y) parameterizes the difference from a generalized Gaussian propagator. The coefficients 
Cj of the low- momentum expansion of l{y) can be related to appropriate dimensionless 
renormalization-group invariant ratios of moments m2j = J^xi^^y^i^)- Let us introduce 
the quantities 
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whose continuum limit is 

t>2j = 1 for a Gaussian critical propagator. One can easily write the coefficients q in terms 
of v*^^: 

C2 = 1 - Vl, 

C3 = l-2vl+ vl 

Ci = l + vl{vl-?,) + 2vl - vl , (30) 

etc. The strong-coupling expansion of G{x) allows to calculate strong-coupling series of V2j. 
Estimates of the coefficients Cj can then be obtained, as we shall see, from the analysis of 
the combinations of V2j corresponding to the r.h.s. of (pOl). 

Another quantity which characterizes the low- momentum behavior of L{y) is the ratio 
s = M'^/Ml where M is the mass-gap of the theory, i.e., the mass determining the long 
distance exponential behavior of G{x). The values s* of s in the critical limit is related to 
the zero yo of closest to the origin: indeed yo = —s*. s* is in general different from one; 
it is one in Gaussian-like models (i.e., when l{y) = 0), such as the large- iV limit of O(A^) a 
models, while no exact results are known at finite A^. 

In the absence of a strict rotation invariance, one may actually define different estimators 
of the mass-gap having the same continuum limit. On the square lattice one may consider 
fi obtained by the long distance behavior of the side wall-wall correlation constructed with 
G{x), or equivalently the solution of the equation G~^{pi = ifi,p2 = 0) = 0. In view of a 
strong-coupling analysis, it is convenient to use another estimator of the mass-gap derived 
from fi: 

= 2 (cosh/i - 1) , (31) 
which has an ordinary strong-coupling expansion 




^' = 75 1 + E«^/51 (32) 



{n has a singular strong-coupling expansion, starting with — In/?). One can easily check that 
Mg/n — i> 1 in the critical limit. Similar quantities and can be defined respectively on 
the triangular and honeycomb lattices, as shown in Apps. ^ and 0. One may then consider 
the dimensionless ratios M^/Mq, M^/Mq and M^/Mq respectively on the square, triangular 
and honeycomb lattices, and evaluate their fixed point limit s*, which by universality must 
be the same for all of them. From the available strong-coupling series of and Mq on 
the square lattice, and Mq on the triangular lattice, and Mq on the honeycomb 
lattice, which are reported, for A^ = 0, 1, 2, in Apps. 0, ^, |C| respectively, we computed the 
ratio M^/M^ up to 16th order, M^/M^ up to 11th order, M^/M^ up to 25th order. For 
the Ising models, using the known exact results for M^, M^, and M^ (see next section), we 
obtained longer series, i.e. M^/Mq up to 20th order, M^ /Mq up to 14th order, M'^/Mq up 
to 29th order. 
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In order to determine s*, and the coefficients q of the low- momentum expansion of L{y), 
we analyzed the strong-coupling series of the ratios M^/M^, M^/M^ and M^/M^, and of the 
combinations of V2j given in Eq. (|30D. Beside the ordinary series in /3, we also considered and 
analyzed the corresponding series in the energy. The change of variable from f3 to the energy 
E is easily performed by inverting the strong-coupling series of the energy E = j3+0{j3^) and 
substituting into the original series in powers of p. We constructed PA's and Dlog-PA's (and 
sometimes as further check also lA's) of both the series in (5 and in E. While PA's provide 
directly the quantity at hand, in a Dlog-PA analysis one gets corresponding approximants by 
reconstructing the original quantity from the PA of its logarithmic derivative. Estimates at 
criticality are then obtained by evaluating the approximants of the /3-series at jSc, and those 
of the i?-series at Ec, i.e., the value of the energy at (3c- In the cases in which Ec is not known 
from independent studies, its estimate may be derived from the first real positive singularity 
detected in the analysis of the strong-coupling series of Xy N = 2, expressed in 

powers of E. 

In our analysis we considered quasi-diagonal [l/m] PA's and Dlog-PA's of the available 
series; more precisely, for a nth order series we considered those with 

11 

l,m>--2, l + m>n-2. (33) 

As a final estimate from each analysis we take the average of the results from the quasi- 
diagonal PA's and Dlog-PA's using all available terms of the series. The errors we display 
are the square root of the variance around the estimate of the results from all non-defective 
PA's indicated by Eq.(|3|). 

By analyzing the above-mentioned series at = 0, 1, 2 we obtained estimates of s*, and 
of some of the coefficients Cj. The results for N = 2, i.e., for the XY model, are reported 
in Table [VI H , those for = 1, in Table [VIII| . Universality among the square, triangular 



and honeycomb lattices is in all cases well verified and gives further support to our final 
estimates. 

For the XY model, the analysis of the £'-series provides the most precise results on all 
lattices considered, leading to the estimates: 

s* = 0.9985(5), 

C2 = -1.5(5) X 10-^ 

C3 = 2(2) X 10"^ (34) 

The errors displayed are a rough estimate of the uncertainty. 

For the Ising model, the two-point function in the scaling region is known analyti- 
cally We obtained a benchmark for our strong-coupling computation by computing 



numerically the two-point function, following Ref. |2^, and performing a numerical integra- 
tion of the results: 



-3 



5 



s* ^ 0.99919633, 
C2 = -0.7936796 x 10 
C3 = 1.095991 X 10~^ 
C4 = -3.12747 X 10~^ 

(35) 
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The analysis of the available strong coupling series on the square, triangular and honeycomb 
lattices lead to the final estimates: 



s* = 0.99908(3), 

C2 = -0.94(4) X 10"^ 

C3 = 1.1(3) X 10"^ (36) 



The agreement with the exact results ( p5|) is satisfactory. But the comparison shows also 
that the errors on the strong coupling estimates of s* and C2, essentially calculated using the 
variance of results from different PA's, are underestimated. We mention an earlier attempt 
to estimate s* for the Ising model by using shorter strong-coupling series on the square and 



triangular lattices 29 



For the self-avoiding random walk model we find: 

s* = 1.0000(2), 
C2 = 0.13(6) X 10-^ 
Icsl < 2 X 10"^ (37) 

The analysis of the coefficients q with i > 3 becomes less and less precise with increasing 
i, but it is consistent with very small values. For instance, we found in all cases {c^l ^ |c3|. 

So, for all considered, our strong-coupling analysis leads to the following pattern of 
the coefficients Cj 

Q < C2 < 1 for i> 3. (38) 

This was also observed in models with > 3 by a study based on large-A^ and strong- 
coupling calculations 0. As a consequence of (|55D, the value of s* should be essentially 
fixed by the term proportional to (p^)^ in the inverse propagator, through the approximate 
relation 

s* - 1 ~ C2. (39) 

This equation is satisfied within the precision of our analysis for A^ = 0, 2, and well verified 
by the exact results of the Ising model, where s* — 1 — C2 — 10"^. 

We can conclude that, like models with A^ > 3, in the critical region the two-point Green's 
function for A^ < 2 is almost Gaussian in a large region around = 0, i.e., \p'^/Mq\ < 1, 
and the small corrections to Gaussian behavior are essentially determined by the (p^)^ term 
in the expansion of the inverse propagator. 

Differences from Gaussian behavior will become important at sufficiently large momenta, 
where G{p) should behave as 

Gip) ~ ^ (40) 
where rj 0: t] = ioi N = 0, and r] = ior N = 1,2. 
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IV. LOW-MOMENTUM BEHAVIOR OF THE ISING MODEL 



So far we considered only the critical limit of the two-point Green's function. It has 
however been known for a long time that the correlation functions of the two-dimensional 
Ising model can be computed exactly for arbitrary values of (3. As a consequence we may 
in principle check directly our computations for every individual coordinate-space Green's 
function. In practice we may perform our checks by exploiting a peculiar feature of the square 
lattice solution: for sufficiently large values of ^Jx"^ + y"^ (in units of the lattice spacing) the 
asymptotic behavior is described by |l30 | 



G{x,y) 



1/4 



{l + z'Y'^l-z') X 



2tT 27r (1 + ^2)2 _ 2z{l - ^2) (cos 01 + cos 02 ) 

where we have introduced the auxiliary variable 

z{iS) = tanh/3. 



(41) 



(42) 



We recognize that the above result (^If) corresponds to the behavior of a nearest-neighbor 
quasi- Gaussian model whose momentum-space propagator has the form 



G{p) 



Z{(3) 



where g{p,P) vanishes at the pole = — M^(/5), 



(43) 



1 - - Az' 



1/4 (1 + ^2)1/2 



and 



2\2 



(l + Z^) 

zil 



4. 



(44) 



(45) 



A straightforward but yet unobserved consequence of this observation is the algebraic rela- 
tionship 



2(coshyUs 



r 



4(coshiyUd 



M\(3), 



(46) 



where /^s and /id are the coefficients of the long-distance exponential decay ( "true mass-gap" ) 
on the side and along the principal diagonal of the square lattice. We verified that Eq. (^6]) 
is satisfied by our determinations of masses from wall-wall correlations and is consistent with 
the known relationship ||2^ 



/is = In coth/? — 2[3. 
We also checked that the residue at the pole fP' = — satisfies Eq. 



(47) 
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Motivated by this piece of evidence we investigated the possibihty that the asymptotic 
behavior of the two-point Green's function of the Ising model on regular two-dimensional 
lattices be always dictated by the structure of the propagator 



G{P) 



Z{(3) 



(48) 



where is the massless (nearest-neighbor) Gaussian inverse propagator appropriate to the 
lattice at hand, and g{p,P) vanishes at the pole = — M^(/3). This conjecture can be 
checked by considering the large- distance behavior of the correlations for the triangular and 
honeycomb lattices, as a function of the direction, and comparing the different available 
mass definitions with each other and with exact results. 

On the triangular lattice we can define a "true mass-gap" fi\ from the asymptotic behavior 
of correlations taken along a straight line of links and a wall-wall inverse correlation length 
/it evaluated in a direction orthogonal to the above defined line (see App. P). In a Gaussian 
model one would obtain 



Ml = - (cosh^/^t 



(^cosh^/^i — (cosh^fii + 2 



(49) 



Starting from the known solution |31 



we checked that the relationship 
expansion of 

M?-- 



2\n i^Vl - z + - y^j - 2 ln(l - ^) - In z, (50) 
is satisfied, since our series for /it reproduces the 



'1 + z^V 1 



1 



(51) 



Finally on the honeycomb lattice two mutually orthogonal inverse correlation lengths 
can be defined by the relationships 



cosh|/iv — 1 
,vl, 



Ml = - l^cosh^/ih - 1 



(52) 



where /iv and /ih are defined from the large-distance exponential behavior respectively of wall- 
wall correlation functions G'^\x) and G^\x) defined in App. 0. The Gaussian relationship 
is 



Ml + 2 



^'Ml + 4Y 



(53) 



Moreover from duality with the weak coupling phase of the triangular lattice model we 
obtained 



73 



v/2cosh2/j-l-l , 

In , ^ h 25 

V2cosh2/:/ -1 + 1 



(54) 
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and we checked that the expansion of the is consistent with 

Ml = ^ (2cosh2/? - 1)^/^ coth/? - 4, (55) 

while Eq. ([53| ) is satisfied to all known orders of the strong coupling expansion. 

In conclusion we may say that the quasi-Gaussian structure of the propagator, described 
by Eq. (^) , is confirmed for all regular lattices and is a remarkable piece of evidence in favor 
of adopting the quantities M^, M^ and Ml respectively as strong-coupling estimators of the 
mass-gap, sharing the property of a well-behaved /3-dependence and of a faster approach to 
universality in models with quasi-Gaussian behavior. It is probably worth observing that, 
since g{0,l3) ^ 0, Eq. (^) does not allow an immediate identification of the moments m2j, 
and in particular M'q ^ M^, and Z{(3) is not the standard zero-momentum wave function 
renormalization but corresponds to the on-shell definition. 
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APPENDIX A: STRONG-COUPLING SERIES ON THE SQUARE LATTICE 

In order to enable the interested readers to perform their own analysis, we present most 
of the series used to derive the results presented in this paper for = 0, 1, 2. This appendix 
is devoted to the square lattice, the following ones to the triangular and honeycomb lattices. 



1. N = Q 



For the self-avoiding random walk on the square lattice, longer series of Mq can be 



obtained from the strong-coupling series of x ^i-nd m2 presented in Ref. pj]. We report our 
series of Mq for sake of completeness. 

= - 4 + 3/? + 2/5=^ + 4/5^ - 10/5^ + 48/?^ - 120 + 368/?^ - 822/?^ + 2008/5^° 

- 4320/3^^ + 10336/3^2 - 22800/3^^ + 56312/3^^ - 129922/?^^ + 327080/?^^ - 768414/3^^ 
+ 1938440/3^^ - 4604254/3^^ + 0(/3^°), (Ala) 

M"^ = 4 + 3/? + 2/3=^ + 4/5^ - 8/5^ + 30/?^ - 52/5^ + 140/5^ - 234/5^ + 596/?^° 

- 1010/3^^ + 2638/?^^ - 4644/?^=^ + 12634/?^^ - 23208/3^^ + 0(/?^^), (Alb) 

V, = ±(3-' + l + Af3+lp^ + 1(3^ _ + lAf37 _ 19^8 ^ 429^9 _ iQ3^10 ^ 511^11 

- 539/3^2 ^ 1468/313 _ 3649/3^^ + ^p^^ - 25668/?^^ + 53|25^i7 _ 154972/5I8 

^ 309|629^19^Q^^20)_ ^^Ic) 
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2. N 



For the Ising model we give strong coupling series which cannot be reproduced using 
known exact results, which are reported in Sec. |IV|. 

/\/f2 _ ^-1 4 I 10 I 134 ^3 I 76 fQ5 , 19394 /q7 q o /o8 , 2070328 /q9 704 ^)10 
IVIq — P - 4 + —p + -^p + YggP + l725"P ~ -^^P + 18711 P 

_l_ 233105490328 /Oil _ 20656 ol2 , 440148292 13 _ 256064 gli , 670306901872438 /ol5 52233344 ,c?16 
~'~ 638512875 P P ^ 729729 ^ 189 ^ ~'~ 162820783125 P 



~'~ 638512875 ^ 45 '^'^ ~'~ 729729 ^ 189 ^ ~'~ 162820783125 P 4725 

I 192016952587260544 ^17 4476104704 ;ol8 , 133522860364557505628 ^19 , r)( ( 

7795859096025 P 93555 P ~'~ 1531329465290625 P ~^ ^\P /' yr^'^d.) 

_ J_/Q-l_|_3,_5_/o,_6^/o3|J9_/i5| 9697 /57 , r, ^^8 _ 339961 /QS , 44 /olO , 8705774291 ^oll 

IqP 4 24/^ 360*^ ~'~ 756-^ ~'~ 37800-^ ~'~ P 37422 3 1277025750 '^'^ 
_ 4514/312 I 3986722469^13 _ 68668 ^14 _ 115832206185781 ^15 , 6752894 ^16 
45 14594580 P 189 P 1302566265000 P ~'~ 4725 P 

21607992820912952 /ol7 , 34100716 /^18 , 71772260149691061407 /ol9 , n( ( ^'?\^\ 

7795859096025 P ~'~ 93555 P ~'~ 6125317861162500 P ^\P )■ y^^U) 



3. N = 2 



p _ /3_|_3^3,1^5_31^7_731^9_ 29239^11 _ 265427 ol3 _ 75180487^15 _ 6506950039 ^17 

-C/ p -\- ^p ^ ^p ^^p ^^qP 1440 y 5040 P 645120 P 26127360 P 

1102473407093 ,019 6986191770643 ,o21 , n( f:i23\ l'A'?fl'l 

2612736000 P 14370048000 P ' /' {-t\Od.) 

X = 1 + 4/? + 12/32 + 34/?3 + + + + '-^p' + ^/J^ ^ 38|93^9 

_l_ 2608499 ,010 _|_ 3834323 ^11 , 1254799^12 , 84375807^13 , 6511729891 ^14 , 66498259799 ^15 

180 '-'^ ~'~ 120 ~'~ 18 ~'~ 560 P ~'~ 20160 P ~'~ 96768 P 

_l_ 1054178743699 |Ol6 , 39863505993331 ^17 , 19830277603399 ,ol8 , 8656980509809027 /ol9 
~'~ 725760 P ' 13063680 P ' 3110400 P ' 653184000 P 
I 2985467351081077 ^20 , 811927408684296587 /^21 , fjfn22\ ('A'^hl 
~'~ 108864000 P ~'~ 14370048000 P ~^ ^\P li \n.OU) 
]\/r2 _ n-l /I , 7/0 , 41/o3 ;o4 , 15,o5 25 , 9491 a7 431 gS , 206411 17803 ^10 

IVIq — P - 4 + 2P + uP ~ P + 16^ ~ YP + '720 P 9~P + 2880 P 360"'^ 

_ 41122019^11 _|_ 876403 ^12 _ 1413373319 ^13 _ 15006841 ol4 , 337093786457 /ol5 _ 4777620367^16 

241920 P ' 1728 P 1935360 P 181440 P ' 130636800 P 1036800 P 

I 17847363647 /ol7 , 68513340691 glS _ 16133717627082721 ,ol9 , r)(Q2Q\ ( A'ir\ 

1741824000 /"^ ~'~ 3732480 P 344881152000 P ~^ ^\P ti \n.O(.) 

= 3-^ -4+^8+ ^8^ - 8* + ^8^ -^8^ + ^8'^ - ^8^ - ^8^ - Mll^io 

ii'-ig ^ ^ -r 2^ ^ ^ ^ 16*^ 6 ^ ^ 720^ 18 ^ 2880*^ 720 ^ 

_ 6120227 ^11 _|_ 24907^12 _ 788579333 ^13 , 95728039 ^14 _ 63069969313 ^15 , r)( f)^^\ ( \'^r\'\ 

241920*^ ~'~ 540 P 5806080 P ~'~ 362880 P 130636800 P ~^ ^\P 1 1 \-t^OU.) 

_ J_/3-l _|_^_|_X/9_|_JL/93 8749 o7 i 67 ^8 _ 122549 o9 _ 2153 ^10 

1^4 i^P 4 32'^ 192*^ 2mP aP 11520*^ ~'~ 48*^ 46080 P 144 P 

_l_ 249335197 oil _ 40951 q12 , 1389732217 13 , 55582271 qU _ 3706449404743 /o 15 , 6252985429 qlG 

~^ 3870720 P 320 P ^ 30965760 P ~'~ 103680 P 2090188800 P ~'~ 2903040 P 

I 75252500337407 /q17 202521546511 /q18 , 163636654204247999 /^19 , /nf /?20N ( k'\p\ 

~^ 27869184000 P 12441600 P ~'~ 5518098432000 P ~^ ^\P >■ \I\0^) 

APPENDIX B: STRONG-COUPLING SERIES ON THE TRIANGULAR 

LATTICE 

The sites of of a finite periodic triangular lattice can be represented in Cartesian coordi- 
nates by 

xih^h) = hVi + km, 
li = 1, ...Li, I2 = 1, ■■■L2, 

r/i = (l,0), ff2=(i^]. (Bl) 
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In order to define a mass-gap estimator, one may consider the following wall- wall correlation 
function 

G^^^(^^l2^=j:G{hv, + kv,). (B2) 

An estimator fit of the mass-gap can be extracted from the long distance behavior of G^\x), 
indeed for x ^ 1 

G^^(x) oce-^*". (B3) 

In view of a strong-coupling analysis, it is convenient to use another estimator of the mass- 
gap derived from fit- 

^ I (cosh^/., - . (B4) 

More details can be found in Ref. 0. 

In the following we show, for = 0, 1, 2, some of the strong-coupling series used in the 
analysis of the O(A^) a models on the triangular lattice presented in this paper. 



1. AT = 



For the self-avoiding random walk on the triangular lattice, longer series of Mq can be 
obtained from the strong-coupling series of x and 1712 presented in Ref. [E3 . 



+ 18356^10 ^ 52792^11 ^ QQ^^Q^U + 63n84^13 + 0{(3'^), (B5a) 

M2 = -4 + fp + 4p' + fp' + fp^ + fp' + ^-fp' + 188P' + + ^P^ 

+ 33181^10 + 0(^11)^ (B5b) 

+ 4|9^10 + 4|3^11 + 21|95^12 + 10(p/jl3 + 0{P'^). (B5c) 



2. N = 1 



For the Ising model we give strong coupling series which cannot be reproduced using 
known exact results. 

A/f2 _ 2/Q-l /I I 32 ;0 I 16 ;o2 I 928 ;o3 , 64 ^4 , 23944 ^5 1648 /o6 , 5008 n7 , 106864 ^8 
IVlfj - 4 -h -g-p -h -g-p -h ^35^ ^ Y ^ ^ ^SS^" ^ 135^ ^ 14175 + ~945~^ 

I 6459424 n9 _ 18680128 alO _ 200433692584 ^11 , 2151999728 al2 , 35136345008 q13 , /n('/?14\ 
~'~ 280665*^ 42525 1915538625 ~'~ 1403325 ^ ~'~ 54729675 ~T~'^\h' ) i 

(B6a) 

- ±/3-l + 3 , 2^ , 1^2 , ^^3 , 4^4 , 2993^5 , 437^6 , ^1^ flT _ 19781^8 
^4-24^ -h 4 -h gP -h gP -h ^35P -h gP -h 567oP + 135P + 14175^^ 945 ^ 

I 403714 n9 I 3986522 /olO _ 354526855073 11 _ 960767707 al2 _ 18090444637 13 , /n('/3l4\ 

280665'^ ~'~ 42525 ^ 3831077250 1403325 54729675 '"'^ >■ l-DUUJ 
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3. AT = 2 



_ 20160371 oil _ 27984359 q12 _ 87289819^13 _ 10256893919 qlA _ 3357272555039 gib , r)( filQ\ 
8640 ^ 4320 5040 ^ 226800 29030400 "-^ J ) 

(B7a) 

X = 1 + 6/5 + 30/?^ + 135/3^ + 570/3^ + 2306/3^ + + + 

_l_ 14339641 q% _|_ 208590287 /olO , 8995595389 oil , 3199713875 /Q 12 , 65793037351 /q13 
30 /-"^ ~'~ 120 ~'~ 1440 "I" 144 840 

I 165647319078571 ^14 , 4600845479023849 ^15 , n( qI&\ ^'R7K^ 

+ 604800 + 4838400 + ^^f^ )' l^""^) 

A/f2 _ 2,o-l /I I «/:?2 I 143/03 I 46/o4 391 /o5 5219 ^)6 5296 7 33287 /pS 

Mn - 3/3 - 4 + + D/3 + + -g-p - -72P - -\^P - 180"^ 

') (B7c) 



604800 4838400 P ^\P )^ 

^2 _ 2,0-1 /I I «/:?2 I 143/03 I 46/o4 391 /o5 5219 ^)6 5296 7 33287 /pS 

^^G - 3P - 4 + + D/3 + + -g-P - -T2P " WP ~ 180"^ 

_ 679729 q9 _ 2052143 qW _ 1436935039 qll _ 4952351659 ql2 _ 87992319949 /Q 13 , ri( all' 
1296 P 1080 P 362880 P 1360800 P 43545600 P , 

i/f2 — 2 q~l _ A \ll q \ «/C?2 I 283 q2, 1 49 /Cj4 _ 53 /o5 _ 8425 q& _ 990757 q7 _ 45549 ^8 
^-'t — 3^ 3 P op 36 P -T 9 P 8 ^ 216 ^ 8640 ^ 160 P 

_ 16833083 q9 _ 35865709 qW , /Of/^llN 
25920 P 25920 ^ "r *^VP ^ 



16833083 q9 _ 35865709 qlQ , /Of/^llN 
25920 ^ 25920 P "i" *^VP ^ 

_ J_/9-l _|_ 3 _|_ 11/0 _|_ 3/o2 _|_ 143/03 I 4/o4 _ i03_/o5 _ .827_/o6 _ 6271/o7 _ 39709 ,08 

^^4 — 24P 4 48P 8P 288P 9P II52P 1728 P 720 P 960 P 

1243813 q9 _|_ 85031 /olO _ 741356239 qll _ 7581779911 ql2 _ 1477616543629 ql^ 
20736 P 3456 P 5806080 P 5443200 P 696729600 P 



(B7d) 

^77616543629 ql^ , n( R^'^\ ^RVo^ 
696729600 P + ^^P )■ 1^"®) 



APPENDIX C: STRONG-COUPLING SERIES ON THE HONEYCOMB LATTICE 

The sites x of a finite periodic honeycomb lattice can be represented in Cartesian coor- 
dinates by 

X = x' + pffp 

X' = hf]i+ l2f]2, 

h = 1, ...Li, I2 = 1, ...L2, p = 0,l, 



(1^), ^-'2=(0,V3), rf,= (l,0). (CI) 



In order to define a mass-gap estimator, one may consider the following wall- wall correlation 
functions 

G^;'\lh)=J2G{iim + i2m), (C2) 

with the sum running over sites of positive parity forming a vertical line; 

G'^\lVSl) - 2^2)^71 + ^2772 +Pf}p): (C3) 

l2,P 

where the sum is performed over all sites having the same coordinate X2- 

Estimators /Xv and /Xh of the mass-gap can be extracted from the long distance behavior 
respectively of (x) and (x) , indeed for x 1 

gI^^x) oc e-'"'^ (C4) 
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In view of a strong-coupling analysis, it is convenient to use the following estimators of the 
mass-gap derived from /i^ and /ih: 

M2 = - (cOSh|/ih - l) , 

Ml = ^ (cosh^/ih - l) . (C5) 

More details can be found in Ref. 0. 

In the following we show, for N = 0,1, 2, some of the strong-coupling series used in the 
analysis of O(A^) a models on the honeycomb lattice presented in this paper. 



1. N = 



For the self-avoiding random walk on the honeycomb lattice, longer series of Mq can be 
obtained from the strong-coupling series of x ^^^^ presented in Ref. [E^ . 



+ 416/3^^ - 1192/3^^ + 2848/3^^ - + 5768/3^^ - ^f3^^ + 20024/3^0 _ 38520/3^^ 

+ 63368/3^2 - 100104/3^3 + 183352/3^^ - Mm^p^^ + 621096/3^6 - 309|i76^27 
+ 1791168/3^^ + 0{(3^^), (C6a) 
Ml = 1/3-1 -4+1(3 + 1(3' + 1(3^ + 1(3'' - 4(3'' + 8(3'' + '^(3" + 8(3'^ + 1(3"' + IQ(3'^ 

+ 1|8/517 + 36^18 + 116^19 ^ l|0^20 ^ 1|2^21 _ 243^22 ^ 4|8^23 _ 3^^24 ^ ^^^25)^ 



75^15 _ 397^16 _|_ 4631 ^17 _ 1051 ^18 _ 349^19 _ 1087 ^20 _|_ 8001 ^21 _ 15883 ^22 _|_ 17989 ^23 
i|91^24 ^ 72^68^25 1 11|63^26 ^ 5Mpi^27 _ HQ^^^f^^S ^ Q^f^29y ' '(C6b) 



2. N = 1 



on _ 1219136 /ol8 
V 945 H 



= -/3~^ _4-|-28/o_124^3i 8576 _ 14692 nl , 5338616 /q9 _ 90947891648 all _ on Q 
IVlQ g/J ^ -t- g /J ^gg/J "t" 2835^ 2025 ~'~ 280665 1915538625 I- 

_l_ 1583805616 /q13 , naO^^ _ 406965884456828 15 _ 608 /ol6 , 360870502928894432 on 
~'~ 7818525 ~'~ 488462349375 ^ 15 ^ ~'~ 116937886440375 

_ 6493740451647884584 ^19 , 624064^20 , 1179228814388026215376 o21 _ 5843265248 o\ 

656284056553125 75 40343570821929375 155925 

_ 48910471162936574893856768 o21, , 405278723648 o2A , 1024764052182397586576416 o2b 

605758715891269565625 "I" 2837835 "r 4754777989727390625 

_ 1697456183968 q2Q _ 950155935558179228231150591072 q27 , 54851554589151328 q2% , /n|'/529^ fPYo^ 
3378375 ^ 1693960980510228821015625 ~'~ 32564156625 ^ '^KH >i K^^'^) 

_ J_/Q-l _|_3_|_X/^ 31_^3, 536 qb _ 3673 ^7 , 667327 /o9 _ 5684243228 |0 11 , o^l2 

^^jJ "T 4 ~'~ 36^ 540*^ ~'~ 2835/-'^ 8100*^ ~'~ 561330'-'^ 1915538625^ "T ^P" 
I 5165551 ql3 _ r-qlA . 93643468635793 ql5 _ 1042 ql6 _ 27120807726815398 ql7 , 643196 ql8 
' 7818525 A^ "I" 1953849397500 15 116937886440375 ^ 945 

I 86312750362752061 ^19 _ 258424 ^20 , 33629885325845121661 ,o21 , 1973399678 ^22 

187509730443750 75 ~'~ 40343570821929375 ~'~ 155925 

_ 6840786318771414403278548 ,o23 _ 106606055168 /o24 , 282271419983204843625526 q25 , 44788447114 ^26 
605758715891269565625 2837835 "T 4754777989727390625 482625 

402797762032926523234583974442 ^27 5992229992104838 ^28 , n( q29\ (C7h\ 
1693960980510228821015625 f 32564156625 y ^ '^Kh' J K^'^) 
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3. AT = 2 



(C8b) 



p _ /Q_i/Q3_|_7/o5_395 /37_|_ 1173 o9 _ 473243 oil , 6293627 /Ol3 _ 346093553 gib 
Hi IJ "t 3/^ 48 "T 40 4320 1^ ' 15120 1^ 215040 1^ 

I 23497364693 oU _ 64962730739719 nl9 , 474090720713083 o2l _ 1641257090013388013 /o23 
~'~ 3732480 2612736000 ' 4790016000 4138573824000 

I 42984420336380838389 a25 _ 11369733294965786406529 ,o27 , 1733398746685522588378351 /o29 , r)( f:i30\ 
26900729856000 1757514350592000 f-^ 65906788147200000 "r J ) 

(C8a) 

^ ^ 1 + + 6/?2 + f /?3 + + ^ 1^6 ^ 1M5^7 ^ 453^8 ^ 6|9^9 ^ 1|^^10 

I 68273 oil _|_ 138307 pl2 , 9157051 ^13 , 42124273 ,ol4 _ 13183321 yol5 , 130286011 alG 

288 360 10080 40320 ^ 645120 ^-"^ 161280 

I 58701184637 ^17 , 246444397309 -0I8 _ 12790078293739 ^19 _ 79551567889 /q20 , 154021837152677 io2l 

8709120 f-^ 43545600 ^ 870912000 f-^ 13608000 ~'~ 1916006400 

I 1452164594591761 io22 _ 393634368786168197 ^^23 _ 4660955848386121 yo24 , 10915691174925870017 io25 
' 28740096000 1379524608000 31352832000 ' 8966909952000 ^ 

I 41989331871750076949 /o26 _ 8481318776641386327367 o27 _ 828979117543657737823 ^28 
~'~ 62768369664000 1757514350592000 329533940736000 

I 5226218120804763962092657 a29 , 21701722199756349611186159 o30 , ^31\ 
263627152588800000 1^ 2109017220710400000 1^ '^XP )j 

A/f2 _ 4/0-1 A I W R 13/o3 I 59/CJ5 r, 3347 iq7 , 35 ^8 , 238009 ^9 493 

IVIq - 3/3 - 4 + -3-/3 - -g-p + - - + -g-/3 + "6480"^ ~ Ts"^ 

_ 19392227 ^11 , 4388 ^12 , 1467214247 /^13 _ 3846767 14 _ 245879581721 ^15 , 362651221 0I6 

181440 45 4354560 12960 195955200 362880 

_l_ 6669774367471 ql7 _ 45385487873 0I8 _ 5401824824719549 ^19 , 68257961593 ^20 

1306368000 10886400 258660864000 3483648 

_l_ 258003704533726433 fj21 _ 362403210060397 fj22 _ 5054778739819764833 ^23 , 72793501494263779 io24 
' 3103930368000 3919104000 15692092416000 ^ ' 172440576000 

I 290687412809274634879279 ,o25 _ 202156253372553206149 io26 _ 43738864245549216552954097 q27 
' 237264437329920000 108637562880000 9490577493196800000 

_l_ 45054828678355702664561 ,o28 , (n(n29\ (CRA 
' 5649153269760000 1 '^iM J; {\^0(.) 

A/r2 _ A a-l A X 10 o 13 o3 , 55 /o5 1^6 7429 /37 2 , 282139 q9 43^10 
^^li, - 3P - 4 + — /? - -g-p + - 3P - - gP + -6480"P ~ 216/' 

_ 26145491 oil I 613 , 2158358071 q13 _ 224587 ,0 14 _ 344839817111 ;ol5 , 1698299 0I& 

181440 P 540P ~'~ 4354560 P 77760 P 195955200 P ~'~ 272160 P 

I 8350838655511 oil , 63590671 /0I8 _ 3061458683224637 /ol9 _ 1894590323 q20 

1306368000 ^ 130636800 129330432000 ^ 52254720 

I 39427276163585267 /j21 , 5154851721889 o22 _ 5303030533425785401 q22. _ 174214610003233 q24 
^ 443418624000 ^ 23514624000 ^ 15692092416000 ^ 147806208000 ^ 

+ O{0^^), (C8d) 

_ J_/Q-l i3,_5_/0 13./33_|__59_/o5|1^6_ 5507 ol _ 71 /08 , 495049 /39 , 961 ^olO 

t^4 12P 4 24P I44P ~'~ I92P ^ 8'^ 4320'^ 96'^ 103680'^ ~'~ 288^^ 

_ 47837987 oil _ 43397 ql2 , 4159541927 ^13 , 11139143^14 _ 643769125241 ^15 _ 1093077841 ^16 

2903040 P 2880 ^ 69672960 <~ ^ 207360 3135283200 P 5806080 ^ 

I 13258559750671 /ol7 , 148926884453 0I8 _ 1148899892904427 ,ol9 _ 6289014713053 -o20 

20901888000 i~ ^ 174182400 P 591224832000 i~ 1393459200 ' 

I 353359211049440273 q2l , 202048769925301 ^22 _ 2704447377391331531 ,o23 _ 278880464307951691 ,o24 
~'~ 49662885888000 P ^ 8957952000 ^ 83691159552000 P 2759049216000 P 

I 54781761414365119518029 ,o25 , 715504863884795060149 q2% _ 10548822433407426077233907 q27 

345111908843520000 ^ 1738201006080000 ^ 13804476353740800000 ^ 

_ 140706622546312581163859 /o28 , 0{ fi29\ (T'Sp'I 
90386452316160000 ^ ^\P >■ 
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FIG. 1. Inx vs. Beside Monte Carlo data from Refs. |14,r^, we show curves constructed 
from the plain series of Inx, from the [10/10] PA of Z^, and from the [6/6/6] lA of biased at 
= 0.5583, such that o ~ 0.50. 
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TABLES 



TABLE L For various values of < 2 and for all lattices considered we report (3c and 7 as 
obtained from a Dlog-PA analysis of the strong-coupling series of aiid (3c^'' and from that of 
^Q. Defective Dlog-PA's, i.e., those with spurious singularities close to the real axis for Re/3 < f3c 
(e.g., Re/3 < l.l/3c) are discarded. An asterisk indicates that most of Dlog-PA's considered are 
defective and the estimate comes just from a few of them, or in the cases where numbers are not 



shown that all Dlog-PA's are defective, so that no estimate can be extracted. 



AT 


lattice 


fl(x) 


7 


Pc 


u 


Y 

4 


triangular 
Eq. # 


0.1728(4) 


1.04(2) 
1.05371... 


* 


* 

0.547925... 


3 
2 


square 
triangular 
honeycomb 
Eq. 


0.258(1) 
0.1875(2) 
*0.339(1) 


0.80(3) 
1.09(1) 
*0.80(5) 
1.08759... 


0.252(1) 

0.1888(4) 

0.3319(1) 


0.65(5) 
0.64(1) 
0.52(1) 
0.574690... 


-1 


square 
triangular 
honeycomb 
Eq. # 


0.3144(1) 
0.2082(1) 
0.42332(5) 


1.13(1) 
1.14(1) 
1.11(1) 
1.15625 


*0.3141 
0.2086(1) 
0.4240(6) 


*0.60 
0.642(3) 
0.64(4) 
0.625 


i 
2 


square 
triangular 
honeycomb 
Eq. d) 


0.34919(6) 

0.2252(2) 

0.48504(3) 


1.233(5) 
1.23(2) 
1.233(3) 
1.23758... 


0.34922(2) 
0.22528(2) 
0.48498(2) 


0.681(1) 
0.687(1) 
0.672(1) 
0.680715... 





square 
triangular 
honeycprnb 
Eq. d) 


0.37900(4) 
0.24087(4) 
0.54117(3) 


1.334(2) 
1.332(5) 
1.341(3) 
1.34375 


0.37905(2) 
0.24092(3) 
0.54116(1) 


0.750(1) 
0.750(2) 
0.748(1) 
0.75 


i 
2 


square 
triangular 
honeycomb 

Eq. d) 


0.40854(1) 
0.25686(5) 
0.59730(2) 


1.494(1) 
1,484(4) 
1.492(1) 
1.49641... 


0.408530(3) 

0.25692(1) 

0.59731(1) 


0.8453(2) 
0.8450(1) 
0.8446(1) 
0.845852... 


1 


square 
triangular 
honeycomb 
Eq. d) 


0.440684(1) 

0.27466(1) 

0.65849(2) 


1.7496(1) 
1.750(2) 
1.750(1) 
1.75 


0.440690(5) 

0.27466(1) 

0.65846(2) 


1.0002(2) 
1.0005(5) 
1.000(1) 
1 


3 

2 


square 
triangular 
honeycomb 
Eq. d) 


0.4804(2) 
0.2967(1) 
0.73371(8) 


2.30(1) 
2.30(2) 
2.313(5) 
2.31987... 


0.4802(2) 
0.2965(1) 
0.7337(3) 


1.31(2) 
1.31(1) 
1.33(1) 
1.33672... 


Y 
4 


square 
triangular 
honeycomb 
Eq. d) 


0.5072(4) 
0.3114(2) 
0.7844(5) 


2.97(5) 
2.91(4) 
3.01(5) 
3.12490... 


0.5066(4) 
0.3111(5) 
0.7845(10) 


1.66(4) 
1.65(9) 
1.74(9) 
1.80413... 


ly 
10 


square 
triangular 
honeycomb 
Eq. d) 


0.529(1) 

0.3230(4) 

0.826(2) 


4.0(2) 
3.7(1) 
4.0(2) 
4.72210... 


* 

0.323(2) 
0.827(4) 


* 

2.2(3) 
2.4(3) 
2.72322... 
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TABLE II. First-order integral-approximant analysis of the 20th order strong-couphng series 
of = (3^^lnx on the square lattice. Asterisks mark defective approximants, i.e., those having 



spurious singularities close to the real axis for Re/3 ^ /3< 



N 


m 


/ 


k 


Pc 


a 


19 


6 


6 


5 


0.5598 


0.55 




6 


5 


6 


* 






7 


5 


5 


0.5563 


0.42 


20 


6 


6 


6 


0.5598 


0.59 




7 


6 


5 


0.5585 


0.51 




7 


5 


6 


0.5565 


0.37 





TABLE III. On the 


square lattice, 


lA analysis 


of the series Ay 


constructed from the 


series of 




and according to CPRM. (Tbiased is obtained by biasing = 1 






N 


m 


/ 


k 






•^biased 


18 


6 


5 


5 


* 




0.522 


19 


5 


6 


6 


1.0082 


0.64 


0.509 




6 


6 


5 


1.0026 


0.55 


0.485 




6 


5 


6 


* 




0.536 


20 


6 


6 


6 


1.0042 


0.57 


* 




6 


7 


5 


1.0043 


0.59 


0.506 




6 


5 


7 


1.0198 


1.05 


0.501 




7 


6 


5 


1.0039 


0.58 


0.523 




7 


5 


6 


1.0035 


0.59 


0.467 



TABLE IV. On the triangular lattice, analysis of the series constructed from the series of 
and according to CPRM, dbiased is obtained by biasing Xc = 1. We noted here that sometime 
biasing the singularity at Xc = 1 gives rise to spurious singularities in the real axis for Xc ^ 1- We 
considered approximants with singularities in the region [0.8,1.2] defective, and they are marked 
by an asterisk. In these cases the estimate of the exponent from non-biased approximants should 



be uiore rclialik;. 



N 


m 


I 


k 




a 


•^biased 


11 


3 


3 


3 


0.9723 


0.192 


0.501 


12 


3 


3 


4 


0.9808 


0.258 


0.504 




3 


4 


3 


0.9881 


0.335 


0.503 




4 


3 


3 


0.9992 


0.490 


0.502 


13 


3 


4 


4 


0.9954 


0.446 


0.510 




4 


3 


4 


1.0017 


0.534 


0.518 




4 


4 


3 


1.0027 


0.549 


* 




5 


3 


3 


1.0026 


0.541 


0.493 


14 


4 


4 


4 


1.0016 


0.532 


0.550 




4 


5 


3 


1.0003 


0.484 


0.477 




4 


3 


5 


1.0015 


0.533 


0.485 




5 


3 


4 


1.0015 


0.533 


* 




5 


4 


3 


1.0021 


0.520 


0.520 




6 


3 


3 


1.0016 


0.523 


0.481 
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TABLE V. The strong-coupling estimates of ^-re compared with some available Monte 
Carlo results on the square lattice, taken from Ref. and obtained by standard Monte Carlo, 
and Ref. by finite-size scaling techniques. The strong-coupling estimates of Cg^ and are 
obtained respectively from [9/9], [10/9], [9,10], [8/11] PA's of Z|, and from [5/6/6], [6/6/5], [6/5/6], 
[5/6/5], [5/5/6], [6/5/5] lA's of biased at (3c = 0.559. We again warn that the errors displayed in 
the strong-coupling estimates are related to the spread of the different approximants considered. 



The asterisk indicates that the number concerns ^^exp) and not ^g- 



/3 








Cg^^' (l5| 


1/2.2 
1/2.08 
1/2.04 

0.5 
1/1.98 
1/1.96 
1/1.94 
1/1.92 
1/1.90 
1/1.88 
1/1.87 


9.320(3) 
18.76(3) 
26.3(1) 
40.3(3) 
52.2(3) 
70.8(7) 

102(2) 

158(3) 

276(11) 

570(30) 

910(60) 


9.318(1) 
18.74(1) 
26.21(2) 
40.08(6) 
52.0(1) 
70.3(2) 

100.7(5) 

156(1) 

269(3) 

549(10) 

860(20) 


9.32(2)* 
18.75(6) 
26.4(2)* 
40.4(4)* 
51.3(9)* 
69.9(8) 


70.4(4) 
100.3(7) 
156(2) 
263(3) 
539(5) 
847(7) 


TABLE VL Summary of the determinations of Pc and a on 
eycomb lattices by different analysis. A bias in the analysis is 
corresponding abbreviation. 


the square, triangular and hon- 
indicated by a subscript in the 


lattice 


series 


analysis 


/9c 




square 


Ix 
Ax 

\^ 
'x 

'f 


DLPA 
lA 

CPRM-DLPA^^^i 
CPRM-IA^^=i 
CPRM-IA^^=i 
PA 

PA 


0.560(2) 
0.558(2) 

0.5579(3) 
0.5583(2) 
0.558(1) 
0.559(1) 


0.53(4) 
0.49(8) 
0.51(4) 
0.50(2) 
0.59(6) 


triangular 


'x 
Ax 

li 

k 


DLPA 
lA 

CPRM-DLPA^^^i 

CPRM-IA^^=i 

CPRM-IA^^=i 

PA 

PA 

IA^ = l/2 


0.3413(3) 
0.33986(4) 

0.3400(2) 
0.3393(1) 
0.3410(5) 


0.52(2) 

0.473(3) 

0.53(1) 

0.50(2) 

0.52(4) 


honeycomb 


'x 
Ix 

k 


DLPA 

lA 

PA 

IA<,=l/2 

PA 

IA^ = l/2 


0.884(1) 
0.877(6) 
0.879(1) 
0.880(1) 
0.878(2) 
0.883(1) 


0.54(1) 
0.4(2) 
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TABLE VII. We report s* , C2 and C3 as obtained from the analysis of the strong-couphng 
series in f3 (first line) and in E (second line) on the square, triangular and honeycomb lattices. 
Beside the spread of estimates from different PA's and Dlog-PA's, the errors displayed take also 
into account the uncertainty on /?c and Ec- We do not report the estimates of C3 by the analysis 



of the /3-series because their uncertainty is much larger than those from the S-series. 



lattice 




s* 


C2 


C3 


square 


13c ~ 0.559 


0.9985(12) 


0.000(2) 






Ec ~ 0.722 


0.9984(7) 


-0.0014(3) 


0.00001(2) 


triangular 


(3c 0.340 


0.9979(11) 


-0.002(2) 






Ec ~ 0.68 


0.9985(11) 


-0.0010(3) 


0.00001(5) 


honeycomb 


I3c ~ 0.880 


0.9988(10) 


-0.001(3) 






Ec ~ 0.77 


0.9987(5) 


-0.0021(4) 


0.00003(2) 



TABLE VIII. For = 1 and A^ = we report s*, C2 and C3 as obtained from the analysis 
of the available strong-coupling series in /3 on the square, triangular and honeycomb lattices. The 
analyses of the corresponding energy series provide consistent but less precise results, so we do not 



report their results. 





lattice 


13c 


s* 


C2 


C3 


N=l 


square 


iln(v^+ 1) = 0.440687... 


0.99909(2) 


-0.00094(3) 


0.000008(5) 




triangular 


\ln3 = 0.274653... 


0.99912(5) 


-0.00098(4) 


0.00001(1) 




honeycomb 


1 ln(2 + V3) = 0.658478... 


0.99907(2) 


-0.00093(3) 


0.000012(2) 


N=0 


square 


0.3790527(2) H 


1.0001(2) 


0.00016(8) 


0.00000(1) 




triangular 


0.240920(1) ^ 


1.0002(4) 


0.0003(5) 


0.00000(3) 




honeycomb 


(2 + V2)-i/2 = 0.541196... §1 


0.9998(2) 


0.00010(7) 


-0.00002(1) 
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